Flavoring the Quark-Gluon Plasma with Charm 
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We demonstrate that a Quark-Gluon Plasma (QGP) with a dilute admixture of heavy quarks 
has, in general, a lower speed of sound than a "pure" QGP without effects from heavy flavors. The 
change in the speed of sound is sensitive to the details of the theory, making the hydrodynamic 
response to "flavoring" a sensitive probe of the underlying microscopic dynamics. We suggest that 
this effect may be measured in ultrarelativistic heavy ion collisions by relating the event-by-event 
number of charm quarks to flow observables such as the average transverse momentum. 



One of the most widely cited findings in ultrarelativis- 
tic heavy ion collisions concerns the discovery of a "per- 
fect liquid" in collisions at the Relativistic Heavy Ion Col- 
lider (RHIC) [IHl. The evidence for these claims comes 
from the successful modeling of the anisotropic expansion 
of the matter in the early stage of the reaction by means 
of ideal hydrodynamics [E-Q, as well as the presence of 
jet-flow correlations that exhibit a conical pattern (ol-fTTj. 

A considerable amount of study has been carried out 
to understand why the transport properties of the system 
created at RHIC are so different (the average 77/s, shear 
viscosity to entropy density ratio is so much lower) than 
the predicted properties of a weakly coupled QGP [I81]. 
This question is currently unsettled, with the observed 
liquid being either described as a strongly interacting 
(t'Hooft coupling constant A ^ 1) QGP |19j. a bound 
state QGP [2^, a (turbulent) Glasma [21| with insta- 
bilities, a thermalizing Hagedorn resonance gas [22j, or 
a "semi-QGP" with Polyakov loops as active degrees of 
freedom |23| . 

Less attention has gone into using hydrodynamics as a 
tool to link QGP phenomenology to observables known 
from first principle quantum chromodynamics (QCD), 
notably the equation of state above deconfinement. It 
should be stressed that studies of collective flow were 
among the earliest predicted observables to probe the 
thermal properties of heated and compressed nuclear 
matter [2J|. As the transverse flow is connected to the 
pressure gradients in the early stage of the reaction, it 
provides information on the equation of state (EoS) and 
might therefore be used to search for abnormal matter 
states and phase transitions [25l - [27j . 

The advent of the Large Hadron Collider (LHC) gives 
us a qualitatively new laboratory for the study of "soft 
flow" observables. A naive extrapolation, based on the 
boost-invariance assumption, logarithmic scaling of mul- 
tiplicity with respect to center of mass energy, and an 
ideal gas equation of state, predicts that the initial tem- 
perature of the LHC will be 1.5-3 times the initial tem- 
perature of RHIC, 500 - 1000 MeV. Thus, one expects 
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that the soft properties of the two regimes are largely 
the same, or at least comparable. However, statistics 
for rare probes produced in hard initial interactions (e.g., 
jets, heavy quarks) should be better by at least an order 
of magnitude and could be used to qualitatively Gauge 
the thermal and statistical properties of the heavy ion 
background. 

In this work we suggest a new way to correlate rare, 
heavy quark probes to the thermodynamic properties of 
the QGP: the change of the equation of state due to the 
presence of a dilute admixture of heavy quarks. We show 
that the speed of sound in a heavy "flavored" QGP de- 
creases with the (small) concentration of charm quarks 
produced in the early stages of the collisions. This leads 
to a novel anti-correlation between the average transverse 
momentum, (pr), and the total number of charm quarks, 
N charm, that could bc mcasurcd at the LHC. 

In general, charm is not expected to be chemically 
equilibrated in heavy ion collisions. The bulk of charm 
content should be produced by "hard" processes in the 
initial state |28l - [3l| at a concentration far above their 
equilibrium expectation [32| . The abundance of cc pairs 
produced in heavy ion collisions at the LHC is expected 
to be reasonably high (^ 10^~^, of course parametrically 
smaller than the total ~ 10* multiplicity), and it follows 
from the strongly coupled nature of the bulk medium that 
they might achieve thermal equilibrium (s^ - lssj . Thus, 
heavy quarks might well function as a "dilute flavor" ca- 
pable of modifying the equation of state of the system 
by a calculable amount, analogously to the way a dilute 
admixture of salt modifles the heat capacity of water. 

Provided the total charm abundance is not changed be- 
tween production and freeze-out (a reasonable assump- 
tion in the dilute limit) and charm is thermally equili- 
brated, the dimensionless quantity p = p/s ~ pT~^ <^ 1, 
where p is the charm number density and s is the en- 
tropy density is a useful parameter to describe the heavy 
flavorness of the event (note that Pc = Pc — yo/2)- If 
the charm distribution is approximately homogeneous, p 
can be assumed to be conserved throughout the hydro- 
dynamic stage provided that charm diffusion is negligible 
(which seems to be the case at RHIC energies 33, 34]). 
Since entropy s is, to a good approximation four times 
the multiplicity [36| , p can be measured event-by-event if 
the detector has appropriate charm reconstruction. Pro- 
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vided an appropriately large sample of events is available, 
the experimentalist can therefore select a sample of ar- 
bitrarily high p "flavorful" events and study the effect of 
flavorness on the equation of state, analogously to the 
way the effect of impurities on a solution's equation of 
state can be studied by adding the required amount of 
impurities by hand in the lab. 

The effects of charm quarks on the QGP EoS have been 
studied within perturbation theory [37], on the lattice 
[ssl . [39| , and also within recently derived sum rules [io*! . 
In the dilute limit relevant for heavy ion collisions, /5 <C 1, 
the effects of charm quarks on the pressure of an inter- 
acting plasma composed of up, down, and strange quarks 
can be calculated to O {T / Mquark) (small at RHIC, and 
most likely at the LHC) by adding a Polyakov loop den- 
sity to the free energy density 



TiT)^MT)+psoiT)FQiT) 



(1) 



where Fq{T) = -Tln£{T) and e{T) is the renormalized 
Polyakov loop (obtainable from lattice calculations) and 
the quantities underscored with denote the values be- 
fore the charm flavor was included (e.g., sq is the lattice 
data for the entropy density of a 2-1-1 QGP [1^). Correc- 
tions to the free energy will generally come with higher 
powers of p and can, thus, be safely assumed to be small 
in the dilute approximation, as well as finite mass ef- 
fects which again are negligible for most if not all of the 
hydrodynamic stage. The density of Polyakov loops con- 
sidered here describes how much the thermodynamics of 
the system changes by the addition of the heavy flavor. 

This simple model predicts that the speed of sound of 
the charmed QGP, Cs{T), is lower than the value found in 
the standard 2-1-1 QGP, Cso{T), near the crossover phase 
transition. This can be seen as follows. The Polyakov 
loop measures the excess in free energy due to the addi- 
tion of an infinitely massive source of fundamental flux 
(see [131 for a discussion of the physical interpretation 
and general properties of Polyakov loops in QCD and in 
holographic Gauge theories). From Eq. ([T]) one can show 
that the entropy density of the charmed QGP is 



s(r) = so(T) 



Uq{T) 



Fq{T) 



(1 



(2) 

where Uq = Fq — TdFq/dT (note that p is tempera- 
ture independent). In confining theories that are also 
asymptotically free, near the deconfinement phase tran- 
sition both Uq and Fq are positive [45|. Thus, the small 
concentration of charm quarks added to the system in- 
troduces correlations that decrease the entropy density 
of the system. The speed of sound is generally related to 
the entropy density s in thermodynamic equilibrium as 
follows 



dlnT 
din s 



Since both Uq and Fq become negative far from the 
phase transition in QCD [S], it is easy to see that the 
specific heat of the charmed QGP increases with p, which 
then implies that the speed of sound Cs{p ^ 0,r) < 
Cs{){T). All this has a simple physical explanation: Cor- 
relations between the medium and slowly moving heavy 
quarks contribute to the energy density but not to pres- 
sure, thus lowering the system's response to perturba- 
tions. Note that its the correlations, rather than the 
quarks, assumed here to be infinitely heavy, that do this. 

Fig. [T] shows our estimate for the speed of sound de- 
rived within the phenomenological model in Eq. ([Ij us- 
ing the expectation value of the Polyakov loop extracted 
from the lattice (2-1-1 QGP with almost physical quark 
masses [i^). One can see that the main effect comes 
from the region near the phase transition (where there 
is a minimum in the speed of sound) but well before the 
Polyakov loop expectation value reaches its asymptotic 
high-T limit, leading to a negative shift of the speed of 
sound from its value in a 2-1-1 QGP (as computed on the 
lattice :4a]). 

An analysis of the recent lattice simulations of a 2+1-1-1 
QGP [1^ validates these conclusions, but the effect we 
are predicting is parametrically larger since it is implic- 
itly assumed in Ref. [39] that p is given by its equilibrium 
value. In heavy ion collisions, initial state interactions 
can bring p up from the equilibrium value by two or- 
ders of magnitude [30l [s^ . |35| . If charm is thermalized 
and dilute, it is reasonable to expect the effect will be 
appropriately amplified. 

Our estimate stops close to Tc, as —T\nl{T) — > oo 
in the confining phase (where the Polyakov loop expec- 
tation value vanishes). Mathematically, one can trust 
our approach as long as the heavy quark is much heavier 
than any other scale in the system, ie — Tln£(T) ^ Mq. 
At some point in the approach to confinement, however, 
this approximation breaks down and the Polyakov loop 
method becomes unreliable. Physically, in a confining 
theory with light quarks, string breaking in the confined 
phase means that the free energy of the system is not 
anymore related to the Polyakov loop expectation value. 

To estimate the contribution of flavoring in the con- 
fined phase, we make the reasonable assumption that, 
just like a flavorless confined QCD thermal system [41], 
flavorful confined QCD is well-described by the hadron 
resonance gas model. In this case, fiavoring can be ap- 
proximated by an admixture of heavy mesons in a gas of 
pions. The latter has a speed of sound of c'/^''* ~ l/v^ 
(ultra- relativistic ideal gas), while the former will have 
a speed of sound of c^™"*' ~ ^JhT / {MImeson) (non- 
rclativistic). It is not difficoult to see that the speed 



of sound of the mixture r42i will go as 
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parametrically smaller than the contribution in the de- 
confined phase, which is just O (p). Thus, the flavoring 
effect on the speed of sound is specific to the deconfined 
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FIG. 1: (color online) The speed of sound dependence on 
the charm quark concentration in the QGP, for p = (cyan, 
dashed lines) and p = 0.1 (red, solid lines). The width of the 
bands denotes lattice uncertainties [131 in the speed of sound 
of a 2+1 QGP, while the thick line denotes the conformal, 
non- interacting value where = 1/3. 

phase, and it is driven by the correlation between the 
quark and the medium due to strong QCD fields, rather 
than just by the large mass of the charm quark. 

The effect of flavoring on deconfined strongly coupled 
QCD is particularly interesting phenomenologically as it 
crucially depends on both the confining and asymptoti- 
cally free nature of QCD. In fact, we expect the medium's 
response to heavy flavorness to be somewhat sensitive to 
the nature of the fixed point present in the ultraviolet 
(Uy). According to the general arguments presented in 
|44l | , in confining gauge theories with a trivial fixed point 
in the UV (such as QCD) M/dT is positive near the tran- 
sition but it changes sign at higher temperatures and, 
thus, the Polyakov loop reaches its asymptotic value at 
high T from above. However, in confining gauge theories 
with a holographic description in terms of supergravity 
M/dT > (this derivative can only vanish at the non- 
trivial fixed point in the UV) . Therefore, the regularized 
Polyakov loop in QCD displays a bump at a given value 
of the temperature in the deconfined phase while in gauge 
theories described via supergravity the loop is a mono- 
tonically function of T that reaches its asymptotic value 
at the fixed point from below Q . 

To illustrate the difference between a theory that is 
non-interacting in the UV and a theory that has a non- 
trivial UV fixed point, we fitted the Polyakov loop data 
[iaj (left panel of Fig. ^ with the function 

expf^-M (5) 

where x = T/Tc, = 0.185 GeV, and ai = 1.24 and 
hi = 2.89. This function, when extrapolated to high T, 
gives an evolution close to what is expected for QCD 
with a peak above Tc We also considered another 

function that is qualitatively consistent with the evolu- 
tion expected for a theory with confinement but with a 



classical supergravity description until T — > oo (i.e, no 
asymptotic freedom) 

/2(x) = exp ^-^^ (6) 

where 02 = 1.75 and 62 = 3.62. Note from the left panel 
of Fig. [5] that both functions fit the data equally well, 
but have a different qualitative behavior when extrapo- 
lated to higher temperatures. As the right panel of Fig. 
[5] shows, the presence of a peak in £ produces a quali- 
tatively different (though quantitatively small) modifica- 
tion of the speed of sound. 

Thus, we expect the medium's response to heavy flavor 
in strongly coupled theories with supergravity duals to 
be, in general, qualitatively different from QCD. We shall 
further elaborate this point. For the A/" = 4 SYM theory 
[451] the modification of the free energy in presence of an 
infinitely heavy quark can be calculated by describing 
the Polyakov loop as a string going from the boundary 
to the black brane horizon. In this picture Fq{T) — 
— ^\/AT. Putting this and the equation of state for the 
■N'sYM = 4 theory [4^1 in Eq. [T]one can show that in this 
approximation = 1/3 regardless of the value of p and 
A, as can be expected for an infinitely massive quark in 
a conformally invariant plasma. 

More realistically, matter fields in the fundamental rep- 
resentation of SU{Nc) can be introduced into the A/" = 4 
SYM theory using D-branes in the gravitational descrip- 
tion [iBi • The effects of heavy quarks on the thermody- 
namics were extensively discussed in '4?] using the D3-D7 
brane setup. The inclusion of heavy quarks in the plasma 
changes the speed of sound and, just like QCD, leads to a 
negative speed of sound dependence on the heavy quark 
content. The change in the speed of sound squared from 

1/3 can be naively expected to be ~ —O {p^j0^- 

The advantage of this approach is that one has full 
control of not just the effects of heavy quarks on the 
equation of state but also on the transport coefficients. 
In this model r]/s = l/Air independently of p [Ssj . 
though both 77 and s get shifted (this is likely to be qual- 
itatively different from physical QCD, for a discussion of 
the dependence of the shear viscosity on the Polyakov 
loop in a weakly coupled framework see p3|). 

These effects produce, at least in principle, observ- 
able consequences that we will now explore. Note that 
p can be easily related to the number of charm quarks 
per unit of rapidity dN charm /dy and the charged multi- 
plicity rapidity density dNch/dy (due to the relationaship 
between entropy and multiplicity, [36|). The first can be 
calculated via perturbative QCD [30| and the second is 
expected to be logarithmic |40|. Thus, using Bjorken's 
formulae [50| 

, _ \ dNcharm/dy 1 NcoU (Tpp~^cc{\/s) 

6 dNch/dy 3 Npart A^l^y In (^^^ 
Using the cross-section shown in [soj, Ncoii ~ ^ti/tt^ ^^^'^ 
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FIG. 2: (color online) Left panel: High temperature extrapolations of the Polyakov loop expectation value, each based on 
functions (see Eqs.((5]and ([6])) that fit lattice data (blue triangles) [43 ] at low temperatures but have a very different high 
temperature behavior, one (black-dashed) compatible with asymptotic freedom, the other expected in a gauge theory described 
by a classical supergravity dual. The corresponding speed of sound, with p = 0.1, are plotted on the right panel. 



the parametrization for multiplicity in [i^ {Aq = 1.14 x 
7r(0.6/TO)^, i?o = 1.41 GeV), we estimate p to be small 
but non- negligible at top LHC energies (see Fig. [S]). 

Of course, this estimate is illustrative only, due to the 
order of magnitude uncertainty in current pQCD calcu- 
lations, as well as the theoretical controversy over dN/dy 
[3l| . Moreover, cc and (to a lesser extent) dN/dy will 
vary event by event. The fluctuation in the number of 
charm quarks is then expected to be Poissonian for a high 
enough event sample, while the fluctuation in dN/dy is 
generally expected to follow KNO scaling [5l|. Hence, p 
is expected to vary considerably event- by- event, a fluc- 
tuation that, as usual, increases for smaller system sizes. 
The crucial issue, though, is that provided that charm 
can be reasonably reconstructed and there is a large 
enough event sample, p is an experimental observable ca- 
pable of serving as a binning class for events (see Fig. 

SD. 

As is well known, there is a connection between the 
speed of sound and the limiting average velocity of a hy- 
drodynamic expansion with shock-like initial conditions. 



(ITVt) freezeout f impart) {Cs) ^ (8) 



where "freezeout" implies averaging over the freeze-out 
hypersurface while the subscript r means the average 
is done over the hydrodynamic evolution. For a shal- 
low shock this result is exact [H^l. While knowledge of 
the initial geometry is needed to establish the form of 
fiNpart), model calculations ,7i, i53i-i55i] indicate that the 
dependence is not washed away even in steeper shocks 
and more complicated initial geometries. 

The final transverse fiow is in return connected to the 




FIG. 3: (color online) p vs ^/s for Pb-Pb collisions computed 
using Eq. (|3]). The dashed lines denote the uncertainties in 
the perturbative QCD calculation [30l ]. 



average transverse momentum 

(Pt) ~T + m {'JtVt) (9) 

In the past (56j . this was proposed as a signature of 
the mixed phase. Since transverse flow, unlike elliptic 
flow, receives approximately equal contributions from all 
stages of the hydrodynamic expansion j57f , the decrease 
of the speed of sound close to Tc (Fig. [1]) could lower (pr) 
for more charmed events with respect to charmless ones. 
The coefficient associated with this heavy fiavoring effect 
would be straightforwardly related to non-perturbative 



Of course, by this we mean "soft" transverse momentum, with a 
~ 1 GeV cut. We underUne this point as to exclude the rather 
high-momentum charm decay products. 
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QCD via Eqs. © and ® , and Fig. □ This effect might 
be easier to measure in smaher systems (such as pp col- 
hsions if they also experience hydrodynamic flow at high 
energies [Hjlllj) due to the greater event-by-event vari- 
ation in p and more precise charm tracking. 

The main requirement of such an analysis is the abil- 
ity to experimentally gauge both the charm quark abun- 
dance and (pt) event-by-event. At RHIC, event-by-event 
charm detection is nearly impossible, since heavy flavored 
particles are reconstructed only from the leptonic decay 
modes ("non-photonic electrons" ), and this branching ra- 
tio only captures 70% of total decays. At the LHC, how- 
ever, it will be possible to find most charm particles in 
each event using primary vertex cuts [goI l6lj | . 

Another obstacle for detecting the admixture's effect 
on the QGP's thermal properties is cross-correlation. 
The average number of charmed particles is positively 
correlated to the average global multiplicity of the event. 
However, (pt) is also correlated to this multiplicity [8j 
due to initial shock depth (the proportionality constant 
in Eq. [8] depends on Npart)- It should be noted, however, 
that the effect of p on (pt) is opposite to the correlation 
with Npart, events with a greater p should have a lower 
speed of sound and hence lower (pr)- Hence, it should be 
possible to disentangle this charming effect from the flow 
response to fluctuations in the initial geometry. Hence, 
and considering that Ncharm ~ NcoiHsions rather than 
N participants, this obstaclc is not unsurmountable pro- 
vided that the event-by-event charm tagging is precise 
enough. In fact, binning tightly in dN/dy and looking for 
a correlation between (pr) and Ncharm should separate 
the (pt) correlation with Npart from the heavy "flavor- 
ing" anti-correlation via the lower speed of sound. 

A possible effect which would give correlations in the 
same sense as the effect proposed here is energy conser- 
vation (roughly, charm quarks need a lot of energy to be 
created, and that lowers (pt))- This effect should how- 
ever be suppressed by factorization and boost-invariance: 
Charm quarks are created from partons with larger 
Bjorken x than the bulk of soft particles at mid-rapidity, 
so the energy they take up comes from regions well away 
from mid-rapidity. If boost-invariant hydrodynamics is 
the correct picture ( diffusion and convection across ra- 
pidity is negligible) we do not expect that energy conser- 
vation will lower (pt) at mid-rapidity. 



While a decrease in (px) with Ncharm could be 
straightforward to measure, it is not the only effect mea- 
surable by binning events on charm content. Since 7]/s 
generally changes with charm content in a Plasma with 
Polyakov loops [l^ f2/e, the ratio of elliptic flow to ini- 
tial eccentricity might depend on the amount of charm in 
the event. Such measurements could be the key to distin- 
guishing between weakly coupled approaches (where t]/ s 
is allowed to vary) and strongly-coupled theories with 
g^ ge-gravity duals (where Ty/s is flxed to leading order 

With very high statistics, this effect could also be seen 
as a variation of the Mach cone angle [9l-[l7| and signal 
height in samples of events where the charm quarks are 
over or under-abundant. 

Finally, hotspots could lead to over-flavored regions 
which would stay together during the hydrodynamic 
stage, leading to a higher transverse momentum fluctu- 
ations (^{Apt)^), and a more anisotropic transverse mo- 
mentum distribution, for more flavorful events. In this 
scenario (pt) when binned by azimuthal angle would be 
correlated to charm direction. 

In conclusion, we have described the modiflcation of 
the thermal properties (the speed of sound and viscos- 
ity) due to the "heavy flavoring" of the plasma by an 
admixture of charm quarks. We have described how this 
effect might be measured in very high energy (LHC and 
higher energies) heavy ion collisions, where a sufficient 
amount of heavy flavor might be produced in the initial 
state to "flavor" the plasma to a level where the change 
in the equation of state will be observable. 
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